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ABSTRACT 

Multi-frequency measurements of the microwave sky can be decomposed into maps of distinct physical 
components such as the cosmic microwave background (CMB) and the Sunyaev-Zel'dovich (SZ) effect. 
Each of the multi-frequency measurements is a convolution of the spectrum on the sky with the bandpass 
of the instrument. Here we analytically calculate the contamination of the component maps that can 
result from errors in our knowledge of the bandpass shape. We find, for example, that for Planck an 
unknown 10% ramp across each band results in a CMB map ST = STqmb — 4.3 x 10~ 3 STsz plus the 
usual statistical noise. The variance of this contaminant is more than a factor of 100 below the noise 
variance at all angular scales and even further below the CMB signal variance. This contamination 
might lead to an error in the velocity of rich clusters inferred from the kinetic SZ effect, however the 
error is negligible, O(50kms _1 ), if the bandpass is known to 10%. Bandpass errors might be important 
for future missions measuring the CMB-SZ correlation. 

Subject headings: cosmic microwave background - cosmology: theory - galaxies: clusters: general - 
large-scale structure of universe 



1. INTRODUCTION 

Small scale anisotropies in the cosmic microwave back- 
ground can arise from a number of sources. In addition 
to the 'primary' anisotropies, generated at the surface of 
last scattering, secondary anisotropies and foregrounds can 
contribute to the observed brightness fluctuations. Sepa- 
rating the components from multi-frequency observations 
is an important part of the data reduction and interpre- 
tation. Here we consider the effect of uncertainties in the 
frequency response of the instrument within the observa- 
tional bands and how that impacts our ability to perform 
■ component separation. 

2. MODEL OF THE SKY 

We assume for simplicity that the intensity in any given 
direction of the sky is the sum of 5 components. The first 
is the CMB itself with specific intensity 
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where B v is a blackbody spectrum and x = hv /krtTcMB — 
z^/56.84GHz is the dimensionless frequency. The k inetic 



Suny aev-Zel'dovich effect (SZ; Sunvaev fc Zc Pdovich [1972 



19801 : for recent reviews see Birkinshaw |1999| and Rephaeli 
1995 ), arising from the motion of ionized gas with respect 
to the rest-frame of the CMB, has the same frequency de- 
pendence as the CMB signal. The second component is 
the thermal SZ effect - one of the primary sources of sec- 
ondary anisotropies in the CMB on small angular scales. 
Ignoring relativistic corrections, the change in the (ther- 
modynamic) temperature of the CMB resulting from scat- 
tering off non-relativistic electrons is 
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where the second expression is valid in the Raylcigh- Jeans 
limit and y is the Comptonization parameter which is pro- 
portional to the integrated electron pressure along the line 
of sight. 

We also include dust, bremsstrahlung (or free- f ree) c mis- 
sion and synchrotron radiation following Knox (199£, see 
below). Specifically for the dust we assume the spectral 
dependence of a modified blackbody with an emissivity in- 
dex of 2 and a temperature of 18K. For the bremsstrahlung 
and synchrotron we assume power-laws in frequenc y wit h 
indices —0.16 and —0.8 respectively (Bennett et al. 1992 ). 
Our model does not include a contribution from extra- 
galactic point sources. For the purposes of this analysis 
we hold the spectral indices of the components fixed (and 
known). Though we do not expect this to be true of real 
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Table 1 

The parameters assumed for Planck. The central 
frequency is quoted in GHz, the beam size in 
arcminutes and the noise (thermodynamic 
temperature fluctuation in a square pixel of side 
'Beam') in fiK assuming a 12 month integration. We 
do not use the 850GHz channel of Planck here since 

OUR FOCUS IS ON COMPONENTS AT LOWER FREQUENCY. 
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astrophysical foregrounds, the impact of this on our anal- 
ysis is negligible. 

3. METHOD 

To understand the effect of an uncertainty in the re- 
sponse of the instrument to a particular signal let us 
postulate the following situation. We imagine that the 
sky is observed at N frequencies, with measurements 9i 
(i = 1, • • • , N). The signal is the sum of M components 
with amplitudes s a (a = 1, • • • , M) such that 



(4) 



where n, is the noise in channel i. We write our linear 
estimator of s a as 



9 Q = W « 



(5) 



Requiring s" to be an unbiased estimator and minimizing 
the rms residual, ((s — s) 2 }, gives 



Pi 



(6) 



where = {riiUj) and F afi = J2ij fiaN.^ 1 fjp is the 
Fisher matrix for the s a . In the case of diagonal noise, 
Nij — erf Si j, the Fisher matrix simplifies to 



F a /3 - ^ fiafip/ erf 

i 

and the estimator can be written as 

ba - <*P a 2 ° l - 

03 i 



(7) 



(8) 



The covariance matrix of the statistical errors for this es- 
timator is ((s a - s a )(sp - sp)} = F~p. 
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Fig. 1.— The angular power spectra, l(l + l)C t / (2-k), of the CMB 
and thermal SZ signals (solid) fluctuations and the CMB and ther- 
mal SZ map noise (dotted). The tSZ contamination of the CMB 
map, CMB(SZ), correlated with itself (with SZ) are shown as the 
dotted curves. We assumed a 10% ramp error on all Planck band- 
passes. 



We are interested in the effects on our component sep- 
aration of a band error which causes the actual bandpass 
to deviate from the design bandpass by 6f. Using Eq. 
with the design bandpass, but replacing 0j with the right- 
hand side of Eq. (Q) evaluated with the actual bandpass, 
and subtracting off the true signal we find 



8s a = (s a - s a ) = Y M c 



P s l3 



where M is the 'component mixing matrix' given by 

M a0 = ^a 7 2 



(9) 



(10) 



and 6fi/3 is the difference between the design and actual 
bandpasses. The bandpass uncertainty then induces a rel- 
ative rms error on component a from component (3 char- 
acterized by 



W a0 
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(11) 



How is 8 fi a related to the bandpass error? If component 
a has frequency dependence g a {v) then 



dv g a (v) 



hi{v) + 8hi(v) 



(12) 



where hi(v) + 8hi(y) is the total bandpass, with the latter 
term the error. It is easy to see that an error in the ampli- 
tude of the bandpass, 8hi{v) oc hi(y), will have no effect on 
component separation since it will 'calibrate out'. In gen- 
eral we shall model this calibration process by demanding 
that 



f i0 = / dv g Q (y) 



hi(v) + 8hi(v) 



1 Vz 



(13) 



where component is the CMB. We normalize g a (v = 
30 GHz) = 1 for all components a. Note that this means 
s a is the amplitude of component a at 30 GHz. 

Note that the above treatment applies in both pixel 
space and spherical harmonic space. In pixel space the 
(Tj should all be calculated for the same pixel size. In 
spherical harmonic space, the Oi (interpreted as errors on 
the beam-deconvolved maps) are ^-dependent: 



<r%{£) = Oi$i exp 



1 



2 V 2.355 



(14) 



where •& and a are the beam and noise defined in Table [l] 
with § converted to radians. 

4. RESULTS 

As an example let us consider simple shifts in the band- 
pass, which we shall model as a top-hat of width Av about 
the central frequency as given in Table [l]. To gain intuition 
let us restrict ourselves to a 2 x 2 subspace consisting only 
of CMB and tSZ signals with a 1GHz shift in the 217GHz 
channel to higher frequency. In this case (at angular scales 
larger than the largest beam) 



M Q0 = 
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-0.00391 
0.00522 



(15) 
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Table 2 

The COMPONENT-MIXING MATRIX, M a p, FOR a 5 COMPONENT model (cmb. sz, dust, synchroton, bremsstrahlung) if all of the 

FREQUENCY CHANNELS HAVE AN UNDETECTED +10% RAMP (SEE TEXT). 
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Table 3 

the normalized mixing matrix, w a g, for a 5 component model (cmb, sz, dust, synchroton, bremsstrahlung) if all of the 

frequency channels have an undetected +10% ramp (see text). 



which indicates that the CMB channel does not 'contam- 
inate' the any of the signals but there is leakage from the 
tSZ channel. That M a o — is a direct result of our cali- 
bration procedure which enforces Sfio = 0. Note that for 
other calibration sources typically used by ground-based 
experiments, e.g. planets, the CMB can, via bandpass 
errors, contaminate other components. 

A more realistic scenario is that the bandpass frequency 
be quite well determined, but the amplitude of the re- 
sponse as a function of frequency be somewhat uncertain. 
This holds for both HEMTs (M. Seiffert, private com- 
munication) and Bolometers (P. Ade, private communi- 
cation). As a simple model of this effect we introduce a 
linear 'ramp' into our otherwise top-hat bandpasses with 
a change in amplitude of 10% across the band. The results 
of including a +10% ramp in all of the channels with our 
5 component model is given in Table |^. 

The elements of the mixing matrix in Table || tell us, for 
example, that our CMB map, Sq will have a contribution 
from dust of — 0.25s2, where S2 is the dust amplitude at 
30 GHz. Fortunately the dust amplitude is very low at 30 
GHz! 

The importance of this contamination is easier to read 
from the normalized mixing matrix, W a p. To calcu- 
late W a p we need to know the rms sky fluctuations of 
the va rious components. We use the model described in 
Knox ( 1999]) and references therein, most notably Bouchet 



& Gispert ( |1999[ ) . The exception is our SZ power spe ctrum 
which we take from White, Hernquist & Springcl ( 2002 ) 
and extend to i < 400 by assuming Cg — C^nn appropriate 
for the Poisson dominated regime (see Fig. ||). 

In Table | we assume s a are the amplitudes of an 
I = 500 spherical harmonic, so that (s 2 ) = Cf. We find 
that all the elements are quite small; the largest contribu- 
tion to the CMB map comes from SZ and the rms of this 
contaminant is .05% of the CMB signal rms. The galac- 
tic contaminants of the CMB map, Wop for (3 = 2, 3, 4 
greatly increase with increasing angular scale. Even so, 



they are very small at all angular scales; at I = 2 Wop 
= -0.006, 0.0002 and 2 x 10~ 5 for dust, synchrotron and 
bremstrahlung respectively. 

In Fig. [l] we show how the SZ contamination of the 
CMB map affects the CMB power spectrum and the CMB- 
SZ cross-correlation power spectrum. Denoting the CMB 
contaminant from SZ as af^ B ^ SZ ^ we have 



C 



CMB(SZ)-SZ 



/ CMB(SZw ; 

( a e m (°. 



SZ 



)*) = M 01 (£)C' ( 



sz-sz 



(16) 

The ^-dependence of the mixing matrix arises from the 
^"-dependence of the noise in the beam-deconvolved maps 
but is quite mild: M (0, 1) monotonically decreases from 
— 4x 10 -3 to its high I asymptote of — 5.3x 10~ 3 . In general 
the cross-correlations are the most affected since they are 
first order in the component mixing matrix. Fortunately, 
for the 10% ramp error, the bandpass error-induced cross- 
correlation is well below the level of the Planck noise. 

Our results are specific to the conservative component 
separation procedure we have assumed. Other methods 
can reduce the statistical noise by including assumptions 
about the statistical pro pertie s of the various components 
(Tegm ark & Efstat hiou ~ 



996; Bouchet, Gispert & Puget 



|l995t Hobson et al. |1998[ ). The results will also differ in 



detail if one adopts a more realistic (and more compli- 
cated) model of the foregrounds. However, we do not ex- 
pect the component mixing matrix to be qualitatively dif- 
ferent for these different procedures or for more realistic 
foreground modeling. The fact that the two-component 
model and the five-component model give similar results 
for both M(0, 1) and M(l, 1) we take as evidence of this 
robustness. 

The bandpass uncertainties can lead to systematic errors 
in the velocities of clusters inferred from the kinetic SZ 
effect. The mixing between thermal SZ and primary CMB 
on a cluster of optical depth r and temperature T e induces 
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an error 

ST ~ -8 x 10- 3 t ( -^U (17) 
\m e c 2 J 

in the CMB signal. If the contamination were all erro- 
neously attributed to kinetic SZ from the moving clus- 
ter, it would bias the inferred velocity by v = 8 x 
10~ 3 (kT e /m e c 2 )c s» 50kms _1 for a rich cluster. This bias 
is negligible compared to oth er sources of uncertainty for 



individual cluster velocities (Haehnelt fc Tegmark 1996 
Nagai et al. 2002. Holder 2002j ) but is comparable to er 



rors that might be achi evable by Planck on b ulk flows in 
10 6 /i~ 3 Mpc 3 volumes ([Aghanim ct al. 200"l|). This sys- 



tematic contaminant of the bulk flows will appear as a 
T e -weighted monopole. Such a pattern is not expected 
cosmologically, would be evidence of bandpass errors and 
could be removed from the data with negligible residuals. 
Note that while we have ignored relativistic corrections to 
the SZ distortion for the purposes of estimating the mag- 
nitude of the effect, they will be impo rtant in the actual 



analysis of the data (Diego et al. 2002) 



5. CONCLUSIONS 

We can conclude that bandpass errors at the 0(10%) 
level are acceptably small for Planck. More sensitive ex- 
periments though may have more stringent requirements 
on the quality of the bandpass measurements because an 
overall scaling of the sensitivity of each channel leaves the 
component mixing matrix unchanged. 
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